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Abstract 

I present a novel class of exactly solvable quantum field theories. They describe 
non-relativistic fermions on even dimensional flat space, coupled to a constant exter- 
7— i \ nal magnetic field and a four point interaction defined with the Groenewold-Moyal 

\ star product. Using Hamiltonian quantization and a suitable regularization, I show 

that these models have a dynamical symmetry corresponding to gl^ © gL^ at the 
special points where the magnetic field B is related to the matrix 9 defining the star 
■ product as B9 = ±1. I construct all eigenvalues and eigenstates of the many-body 

Hamiltonian at these special points. I argue that this solution cannot be obtained by 
, any mean-field theory, i.e. the models describe correlated fermions. I also mention 

> \ other possible interpretations of these models in solid state physics. 

1. Introduction. Exactly solvable models play a special role in quantum field theory. 
They provide a testing ground for general methods and concepts applicable to other, 
^ ■ more realistic, models which can be studied only in approximations. Moreover, such 
approximations usually are based on some exact solution. The majority of the known 
exactly solvable quantum field theories are free (non-interacting), and most of the few 
known interacting ones are in low dimensions and/or have supersymmetry (see e.g. Jl|, [|). 

In this paper I find and study a particular class of interacting field theories without 
supersymmetry (as far as I can see) and which are exactly solvable in 2n + 1 dimensional 
spacetime with n = 1, 2, . . . arbitrary. They describe non-relativistic fermions coupled to a 
constant, external magnetic field and four-point interactions defined with the Groenewold- 
Moyal star product .[] These models are specified by two antisymmetric 2n x 2n matrices 
B and 9 which define the magnetic field and the star product, respectively (precise defini- 
tions will be given below). Using the Hamiltonian framework, I will solve these models for 

1 To avoid misunderstanding I stress that the free part of the model is the standard Hamiltonian of 
fermions in an external magnetic fields: No star product is used there. 
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the special cases B9 = ±1 and arbitrary interaction constant. A crucial point in the argu- 
ment is to expand the model in a convenient basis which brings it to a matrix form which 
makes manifest a huge dynamical symmetry. It is interesting to note that, in this latter 
form, this model also has a natural solid-state physics interpretations: it also describes 
correlated fermions in 2 + 1 dimensions M . Readers interested only in this aspect of this 
work may skip all the discussion on non-commutative spaces and start immediately with 
Section 4 below. Another aspect is the mathematical structure involved in the solution: 
The above-mentioned dynamical symmetry corresponds to the Lie algebra gl^Qgl^, and, 
from a mathematical point of view, the solution boils down to decomposing the regular 
representations of all permutation groups into irreps. The latter is a classical problem 
in group theory, and I can solve the model by using some beautiful group theory results 
obtained a long time ago (books containing proofs of the results needed are, e.g. [|, |[| [7]). 
However, I will first derive the solution by an explicit construction exploiting the dynami- 
cal symmetry. That method is perhaps less elegant but is, as I hope, a good starting point 
for getting an intuitive understanding of the solution. (The former elegant method is, to 
my opinion, only easier for experts in group theory. One could also read the two methods 
together as a pedestrian proof of certain group theory results.) The group theory point 
of view suggests further generalizations of the model with additional (peculiar) p-body 
interactions where p = 3,4, . . .. These additional interactions correspond to higher order 
Casimir operators and leave the model solvable. 

As mentioned, the models considered can be naturally interpreted as field theories of 
fermions on a noncommutative phase space (time remains commutative): using the star 
product implies that the components of the spatial variables x = (x 1 , . . . ,x 2n ) are made 
noncommutative 

X" * a? - a? * = -216!™ (1) 

{n,v = l,...,2n). Thus, generalizing a standard field theory model by replacing the 
pointwise product of fields by the star product can be interpreted as making space non- 
commutative. In a similar manner, coupling the fermions to a magnetic field can be 
interpreted as a further generalization where momentum space is made noncommutative 
as well: one replaces the commutative momentum operators by noncommutative ones, 
— z<9 M — > — id^ — B fJLV x u . Field theory models with noncommutative space have received 
much attention recently (for review and a fairly exhaustive list of references I refer to 
H ^ [TtJ). In particular, these standard noncommutative field theories are known to 



simplify, become soluble, for the limiting case 9 = oo [JTT], [L2| . I show here that there are 
other solvable cases in the extended family of field theories with noncommutative phase 
spaces, namely all those where B9 = ±1. It is interesting to note that 9 = oo is precisely 
the limit where the kinetic energy term is irrelevant, whereas the solutions obtained in 
this paper allow to study the effect of (a particular kind of) kinetic energy as well. 

It is worth noting that the field theories considered here are special in that they look 
the same in position and in Fourier space and, in particular, for B9 = ±1 they are 
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(essentially) invariant under Fourier transformation [13]. I believe that it is the latter 



property which makes this model exactly soluble: the boson models discussed in [|13|] 



seem to be soluble at these special points as well |TJ]]. I also mention that there is a 
(classical) field theory model of bosons which somewhat similar to the one we solve and 
whose integrability was noted already in jnj. I believe that the quantum version of this 



model should be solvable by similar methods as the ones used here. 

It is interesting to note that these field theories in 2 + 1 dimensions may provide 
interesting models for a quantum Hall system: I add to the standard free quantum Hall 
Hamiltonian a particular four-point interaction of the fermions, and the resulting model 
is exactly solvable. In fact, the model I solve is somewhat more general in that there is 
also a term describing a confining electrical potential. I will elaborate this interpretation 
slightly in the conclusions, but I should stress already here that my results (i.e. eigenvalues 
and eigenfunctions) provide only a first step to understand the physics of this model. 

The plan of the rest of this paper is as follows. In the next section I first define the 
models in the context of noncommutative field theories (Section 2). In Sections 3-6 I solve 
these models in 2+1 dimension. The key is two mathematical facts about Landau eigen- 
functions which allow to write the Hamiltonian in a matrix form (Section 3). In Section 4 
I give a precise meaning to the model using a natural regularization and constructing its 
Hilbert space representation. In this Section I also demonstrate the dynamical symmetry 
of the model. I then give two methods of solution, one using a somewhat pedestrian 
approach (Section 5), and another using group theory (Section 6). The generalization of 
this to 2n + 1 dimensions in Section 7 is remarkably simple: from an abstract point of 
view, changing the dimension does not make much of a difference. Section 8 contains a 
short discussion of generalized models with additional p-body interactions (p = 3,4. . .). 
I end in Section 9 with a few remarks on possible applications, generalizations, and open 
questions. For the convenience of the reader two appendices are added: Appendix A con- 
tains an elementary proof of the above-mentioned facts, and Appendix B is to exemplify 
the general group theory results. 



2. Definition of the models. The models we consider are defined by a Hamiltonian 
H = H + H- m t where the free part is H = J d 2n x^(x)H B ^(x.) with 

H B = {-id, - B, u x v f (2) 

the (generalized one-particle) Landau Hamiltonian; B = (B^) is the matrix defining an 
external magnetic field. (Our conventions for B and 9 differ from the usual ones by factors 
of 2, in order to simplify some formulas. In particular, the magnetic field is 2B.) The 
interaction is 

H int = go J d 2n x (*t * * * *t * * ) ( x ) (3) 

where * is the Groenewold-Moyal star product as usual (see e.g. Eq. (9) in fl3]). The real 
coupling constant g is arbitrary. 
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Formally, I define the quantum field theory by postulating standard anticommutator 
relations for the fermion fields: 



tt(x)*t(y) + ^(y)^(x) = 5 2 "(x - y) 



(4) 



etc. I will give a precise meaning to this model by introducing a particular cutoff A taking 
care of all potential divergences. The limit A — > oo turns out to be rather trivial, at least 
if one is only interested in the eigenfunctions and eigenvalues of the Hamiltonian. In fact, 
we will be able to solve a more general class of models with 



where we allow for an arbitrary linear combination of the Landau terms with B and —B 
(a,b > are arbitrary constants), and I also find it convenient to inserted a chemical 
potential fi (= arbitrary real constant). This generalization is interesting since, as is well- 
known, the Landau Hamiltonian Hb is highly degenerate, but in this extended family 
of models we can also study the case where this degeneracy is lifted. It is worth noting 
that, in the context of a quantum Hall system, the second term has a natural physical 
interpretation as a confining electrical potential. As we will see, there is only a single 
divergence in the model which can be removed by normal ordering or, equivalently, an 
additive renormalization of /i. 

We now rewrite our model by expanding the fermion fields in a convenient basis of 
one particle wave functions. The resulting Hamiltonian has a structure such that we can 
construct a complete set of exact eigenstates and compute explicitly its partition function. 
To simplify the presentation we first concentrate on the special case In = 2 and then give 
the generalization to arbitrary dimensions 2n. 

Remarks: To prepare for the results in the next Section it might be helpful to point 
out a complimentary interpretation of the 'noncommutative phase-space' model above: 
As known since a long time, the star product corresponds to a representation of Hilbert 
space operators by functions ||. Thus the non- commutative field theory amounts to 
replacing the one-particle states of the theory by Hilbert space operators: plane waves 
are replaced by operator obeying 



where k-^k' = k^O^k^. The latter relation is equivalent to interpreting the components of 
x as operators obeying the commutator relations [x^, x u ] = —2i6 fiu (the Fourier variables 
k' v are real numbers). Thus, instead of fields ^(x) multiplied with the star product, 
we could also use fields 




(5) 



i(k+k')-Xgik-0k' 



(6) 




(7) 
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with 'standard' products (and the Fourier transform \1/ an ordinary function). Indeed, for 
a = 1 and 6 = one could also write the Hamiltonian above as 



U = Trace i2(R2 ,) (^( x ) p 2^( x ) + g[& (±)} 2 ) (8) 

with 



x 



af-ie^dv, Pr = -id li -B llv 3? (9) 



the 'Schrodinger representation' of the non-commutative field theory. Eq. (|SJ) implies 

^,f] = -1(^ + 5^), (10) 

i.e. if B9 = —I the 'non-commutative positions' x commute with the 'non-commutative 
Laplacian' p 2 ! This suggests that the model should be special at B9 = —I, and in 
particular it should have a huge gauge-like symmetry. 

We also recall that, in a standard field theory, one can expand the fields $ in a 
basis \£) (Dirac-notation) of eigenfunctions parametrized by positive integers £, and the 
expansion coefficients are infinite vectors (Ai)^ 1 . The corresponding basis \£){m\ for the 
Hilbert space operators is labeled by two integers £, m, and we thus should expect that 
fields can be represented by infinite matrices (Ae m )f m=1 . This suggests a close relation 
of noncommutative field theories and matrix models (various such relations have been 
discussed before || §, 10 1, but I believe the one used in the next Section is particularly 



simple) . 

These remarks will be made precise in the next Section. 



3. Matrix form of the 2 + 1 dimensional model. We assume 2n = 2. My discussion 
will be based on two mathematical facts. Both facts are known since many years in the 
context of phase space quantization (a recent discussion including references to early work 
can be found in |IB|). More recently they have been used in the context of noncommutative 
solitons |TT| (for other references see || |9], [TtJ). However, I would like to stress their 



central importance for noncommutative field theory somewhat more than usually done in 
this context (not only for the models discussed here but in general). I have formulated 
them such that they are true, as they stand, also in 2n dimensions (this is shown further 
below). For the convenience of the reader, (elementary) proofs of these facts are given in 
Appendix A. 

Fact 1: There is a complete, orthonormal basis of one-particle wave functions <^ m (x) ; 
labeled by positive integers i and m, and which have the following star product relations, 

4>lm * <t>i'm' — r <5m,£'</W (11) 

with some positive constant r.^ Moreover, <p\ m = (p mi . 



2 This is true in any dimension In. For completeness I quote r 2 — (47r)™-\/dct((9). 
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This suggests to expand the fermion fields in this basis, 

v&(x) = ^^ m 0, m (x), *t (x ) = ^ i4 t n ^t m ( x ) j (12) 

where the fermion operators A^ obey the usual anticommutator relations, 

^im^i'm' + ^I'm'^lm, = °~l,£ ,0 ~m,m' 

A simple computation (using Fact 1 and J d 2 x 4>\ m -k <j)£> m i (x) = <5^'5 m , TO ') yields, 

= E^i^^^ft = ^Tr^AAU) (14) 

where g = go 7 " 2 - Here we interpreted the A lm as components of a (infinite) matrix A 
with adjoint A^ defined as (A^)e m = A^ ml , and Tr is the usual matrix trace (sum of the 
diagonal). Thus this basis (p£ m allows us to write the interaction in a matrix form. We 
now observe that this basis has a remarkable physical interpretation. 

Fact 2: The functions 4>i m are common eigenfunctions of the Landau Hamiltonians H B 
and H_b in (Hi for B = O^ 1 . The corresponding eigenvalues Eg and E m only depend on 
£ and m, respectively. 

The latter property is the well-known degeneracy of the Landau Hamiltonian. We also 
recall that the eigenvalues are identical to those of a harmonic oscillator, 

Ei — 4\B\(l — (15) 

(recall that we label states by I = 1,2,...) and similarly for E m . Thus, if we choose the 
magnetic field as B = the free part of the Hamiltonian in @ has the following simple 
form, 

H, = Y,{E m + E,)A\ m A lm (16) 

m,i 

where we set a = 1 and defined Eg = bEp — fi. The latter is only to simplify notation: in 
our solution below the explicit form of E m and E(, is not needed. 

Remarks: Fact 1 shows that the functions rtp^ m provide a representation of rank-one 
Hilbert space operators \t)(m\ (Dirac bra-ket notation). The existence of such functions 
should not be surprising (see the remarks at the end of the last Section). What I find 
remarkable, however, is that these functions are old friends to anybody familiar with the 
theory of the fractional quantum Hall effect (Fact 2) WJ\ . 

4. Regularization and integrability. We consider the eigenstates and eigenvalues of 
the field theory Hamiltonian 7i = TCq + Hint defined in Eqs. fll6|) and ([14]) , with fermion 
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operators obeying the usual anticommutator relations (13). We use the regularization 
defined by restricting the quantum numbers £, m — 1, 2, . . . A with A < oo. It is interesting 
to note that A provides a UV (= short- distance) and IR (= long-distance) cutoff at the 
same time (the interested reader can find a more detailed discussion on this in JIB]). 
The regularized model is defined on the fermion Fock space J-\ generated by the fermion 
creation operators A\ m , 1 < £, m < A, from a normalized vacuum Q defined by 

A £m tt = for all ^m, (17) 

and such that f is the Hilbert space adjoint. Since T\ is finite dimensional^ all potential 
divergences are taken care of by this regularization. 

It is convenient to introduce normal ordering : • • • : in the interaction as usual.Q We 
observe that 

:H int : = n- mt -gAj2 A li A mi, (18) 

i.e., normal ordering amounts to a shift of the chemical potential, /i — > — gA. This shift 
diverges as A — > oo and corresponds to a renormalization of the chemical potential. After 
this renormalization all eigenvalues and eigenstates have a well-defined limit A — > oo. To 
see that we recall the following natural basis in the fermion Fock space JFa, 

\N) = A\ im A\ 2m2 ...A, Nm n (19) 

distinguished by the fermion number N = 0,1,2,... and labeled by N distinct pairs 
(£j,rrij). As we make explicit below, the Hamiltonian 7i = Ho+ '■ 7~L m t : acting on all 
\N) is always well-defined even in the limit A — > oo. Thus normal ordering is enough to 
remove all divergences in all the eigenfunctions and eigenvalues of the model. 



In the following I write Ti, m t short for : 7i 



hit 



Remark: The computation of other quantities of physical interest may still require some 
regularization, of course. For example, interpreting the model as a quantum Hall system, 
one would be interested to study the model at fixed fermion density, i.e compute the 
partition function of the model at finite cut-off A, and then take the limit A — >• oo at fixed 
expectation value of N/A. 



Before constructing the solution we give an indirect argument showing integrability of 
the model. For that we define the operators 

Plm ~ A kt A kmi Plm = A \k A mk (20) 



3 Tk = with TV = 1 b 



4 • A A^ A ■ — A^ A^ A A 

■ Ii m 1 e 1 ^ l m 1 i2^ ± m 2 £2 rn 2 li ' — ^mjij ^m^i tnifc^i ' 
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and observe that they provide two commuting representations of the Lie algebra gl^: using 
the fermion anticommutator relations one can show that the p's obey the commutator 
relations 

\Plmi Pt'm') — bmfiPtm 1 ~ ^m'/Pi'm (21) 

and similarly for the p's, and [pe m ,Pe'm'] — (the latter fact is not obvious but follows 
from a straightforward computation). Moreover, 

Pin = Pme ( 22 ) 

and similarly for the p's. We thus see that the operators p and p represent the Lie algebra 
©gU- We also note 

pirrSl = ptmft — for all £, m. (23) 



Remark: For the readers appreciating fine points in analysis I should stress that when 
I write gL^ I mean the inductive limit A — > oo of gl A with A the matrix cutoff described 
above: since we only consider the action of operators on finite particle vectors \N), this 
limit is trivial. Put differently: it is true that some of our Fock space operators become 
unbounded for A — > oo, but we always consider them on a particularly nice common 
dense invariant dense domain. (One could without difficulty replace everywhere in our 
discussion 'oo' by 'A'). 



I now discuss the dynamical symmetry of the model. Note that the free part of 
the model is a linear superposition of Cartan elements of these representations, and the 
interaction is proportional to a Casimir operator. More specifically, (0) is equivalent to 

7~(-0 ^ ^ j y^mPmm ^mPmm) i (24) 
rn 

and the interaction can be written in the following two equivalent forms, 

T^int = 9 2^ : PimPmi '■ = ~9 : PlmPml '■ (25) 

(the first equality here is obvious from (fl^) , and the second is obtained by interchang- 
ing the two A's in ( |HD using that all fermion operators anticommute under the normal 
ordering symbol). This shows that Ti is a sum of commuting operators. Moreover, even 
though the p's and p's commute with the interaction ?^i nt , they do not commute with 7i : 

[Ho :PU'\ = {Eg, — Ef ) pw 
[HoiPmm'] — {E m — E m i)p mm i . (26) 
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This shows that the Lie algebra gl^ © gl^ is a dynamical symmetry for our model. This 
rich symmetry structure suggests that group theory should provide powerful tools to 
elegantly solve this model. We also note the following commutator relations 

[P£>i»'A em ] = h",iA\, m 

[P m ' m "'\m\ = bm",m,A\ m , (27) 

which will be useful for us below. 

We now compute the action of the Hamiltonian on the vectors \N) defined in ([19]). 
All \N) obviously are eigenstates of the free Hamiltonian 7{q in (|T6| ) with eigenvalue 

N 

S = Y,(Ei J +E m] ). (28) 

3=1 

Moreover, by a straightforward computation using the fermion anticommutator relation 
we obtain the following equation 

H- mt \N) = 2g T w\ N ) (29) 

l<j<k<N 

with Tu k \ a transposition operator defined as follows, 

T(jk)A eimi ■ ■ ■ A iNmN Q = A iimi ■ ■ ■ A £ , mk ■ ■ ■ A lkm . ■ ■ ■ A iNmN tt (j < k) (30) 

(i.e. rrij and rrik are interchanged). Eqs. (|29|)-(|30D will be the key to our solution. Note 
that the operators T( jfc ) generate a (highly reducible) representation T of the permutation 
group SV which acts on the states \N) by permuting the quantum numbers rrij. Further 
below we will show how to use the representation theory of Sn to construct eigenstates 
and eigenvalues of the model. However, we first turn to a different approach exploiting 
the dynamical symmetry. 

Remark: A physical interpretation of the relations (|26|)-(p7D is as follows. The one particle 
energies of our model are sums of two parts, Ei m = Eg + E m . Since \Hq, Al m ] = (Eg + 
E m )A tm , applying the fermion operator A\ m adds the corresponding one-particle energy 
to the state, as usual. The peculiar feature of our model is that we have operators 
and p mm i allowing to change only parts of the one-particle energy. 

5. Solution I. Pedestrian approach. We observe that we can generate new eigenstates 
of Ti from known ones by applying operators pw and p mm , r - according to our discussion 
above this gives new eigenstates where the eigenvalues are changed by amounts E m — E m i 
and Eg — E^, respectively (this follows from (BS)). We thus can obtain many eigenstates 
from a sufficiently large number of special ones. 
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We first construct special states of the form (|T9|) which are eigenstates of all transpo- 
sitions Tuty and thus are trivially also eigenstates of 7ii nt . Obvious such states are those 
where all rrij are the same, e.g. 

^1,1^2,1 ' ' ' ^lv,i^ 

(note that due to the Pauli principle, i.e. (A\ m ) 2 = 0, such a state is non-zero only if all 
£j are different): since all rrij = m k , applying Tufy does not change anything, i.e. this is 
an eigenstate of all T^) with eigenvalue equal to 1. Thus this state is eigenstate of Hj nt 
with eigenvalue 2gc where 

c = (number of T (jk) with 1 < j < k < N) = ^N(N - 1). 

A more general such state is 

| [Ax, A 2 ]> = AlA,! ■ ■ ■ KAA,2 ■ ■ ■ A l,2^ 

with Ai + A2 = N, i.e. we have two groups of fermion creation operators where the 
rrij are equal within each group but different in the other group. This is obviously an 
eigenstate of all Ttjk) with j and k both in the same groups (i.e. 1 < j < k < Ai and 
Ai + 1 < j < k < Ai + A2), and the eigenvalue of all these is +1. Moreover, it is also in an 
eigenstate of T 1Al+1 but with eigenvalue equal to —1: interchanges A\ x and A\ 2 , 

but due to the fermion anticommutator relations this is the same as multiplying with — 1. 
Similarly, all I},a x +.j have eigenvalues —1. However, applying Ti } \ 1+2 changes A\ x to A\ 2 
and Al 2 to A\ v Now at least one of the operators A\ 2 or A\ x appears twice, and the 
resulting state thus is zero (Pauli principle). The same is true for all other Tu^\ with j 
in the first and k in the second group and £ k 7^ if all these annihilate the state |[Ai, A2]) 
due to the Pauli principle. We thus conclude that |[Ai, A2]) is an eigenstate of Ttmt with 
eigenvalue 2gc where 

c = -Ai(Ai - 1) + ^A 2 (A 2 - 1) - min(A 1 , A 2 ) 

(we counted the number of T(jk) belonging to the same group and subtracted the number 
of T(j fc ) belonging to different groups but with £ k = £j). It now is obvious how to generalize 
this to states with an arbitrary number L — 1,2, ... N different groups: The crucial point 
in the previous example was not only that the rrij in each group coincide, but also that a 
maximal number of the £j in the different groups are the same. We thus define 

I [A]) = A 1)1 A 2)1 ■ ■ ■ A Xi l A 12 A 2 2 ■ ■ ■ A\ 2 ^ 2 ■ ■ ■ A 1L A 2 L ■ ■ ■ A\ l L Vl (31) 

with a set [A] = [Ai, A2, . . . , Ai] of positive integers A, such that Ai + A2 + . . . + Xl = N. 
Similarly as above we check that this is an eigenstate of all Tj k with eigenvalue +1 if both 
j, k belong to the same group (i.e., rrij = m k ), —1 if j, k belong to different groups but 
£j = £ k , and in all other cases. Thus |[A]) is an eigenstate of 7i int with eigenvalue 2gc[\] 
where 

L 1 

C W = Yl 2 A ^ Al ~ ^ ~~ min ( A ^ A i) 

1=1 l<i<j<L 
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We thus have obtained various particular eigenstates of Ti with different eigenvalues of 
the interaction H.- m %. To avoid confusion we note that each of these states represents a 
whole class of states, e.g. using ( p7| ) we can write any eigenstate of 7i where all rrij are 
the same as follows, 

A im A l 2 m ■ ■ ■ A \ N nP< = Pi^PliZ ■ ■ ■ P£ N APm,l) N \ W}) , 

where \[N]) = A\ x A\ 1 ■ ■ -A' N1 Q is the particular state which for now is the only one of 
this kind taken into account. In fact, at this stage we can further restrict ourselves to 
states | [A]) which are in one-to-one correspondence to partitions [A] = [Ai, A2, • • • , \n] of 

N, 

Ai > A 2 > ... > A L > 0, X 1 + X 2 + --- + X L = N. (32) 
This allows us to write a somewhat simpler formula for the possible eigenvalues of 7i int , 

L j L 

c [A] = ^2 Aj(A4 + 1) "S zA - (33) 

i=l i=l 

All these states |[A]) are obviously also eigenstates of TCo with eigenvalue £® = £q as 
in (|28|) but with £j = £j and rrij = rhj where 

(44 = (l,2,...,A 1 ,l,2,...,A 2 ,...,l,2,...,A L ) 

(mi,m 2 , ■ ■ .,m N ) = ( 1, 1, . , 1 , 2,2,^. . , 2 , . . . , L,L, . ..,L) . (34) 

Ai A2 Xl 

Using this notation we can also write 

I [A]) = A\ _ A\ _ ■■■A\ _ n (35) 

Thus I [A]) is an eigenstate of TC with eigenvalue £° = £q + 2gc[x\. As discussed above, we 
can now generate other eigenstates with different eigenvalues by applying operators 
and p mm i to known eigenstates. In particular, the state 

^ = Pe 1 IiPe 2 /2 ' ' ' Pe N ,£ N Pmi,miPm2,fh2 ' ' ' Pm N ,fh N \ [A]} (36) 

is an eigenstate of H but with eigenvalue 

£ = £ + 2gc [X ], £ in (28) and c [X ] in (||) (37) 

(to see this, apply Ho to $ in (|36|) , move it to the right using repeatedly the relations 
in (|26"1), and note that the energy differences thus generated add up to £ — £®.) It is 
worth noting that the second of the quantum numbers for the operators p and p in (|36|) 
were forced on us by Eqs. (p7|) and (p3|): any other choice would have given zero (this 
will become obvious further below). I will argue below that we thus have obtained all 
eigenstates and eigenvalues of our model. In fact, we obtained too many: many of the 
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states in fl3~6|) are actually actually zero and thus not all of the eigenvalues in fl3~7|) exist. 
For example, if all rrij = m are identical, $ in (^) is nonzero only for [A] = [N], i.e. in this 
case only the eigenvalue £ fl57|) exists where c = iV(iV — l)/2. To characterize all non-zero 
states is a non-trivial combinatorial problem which we will discuss in more detail further 
below. In the rest of this Section we derive a more explicit formula of the eigenstates $ 
(0)- 

It is useful to first consider the special case [A] = [N], i.e. 

$ = Pli,lPti,2 ' ' ' Pe N ,Npm 1 ,lPm 2 ,l ' ' ' Pm N , 1-^1,1^2,1 ' ' ' ^,1^' 

We commute p mjVi i to the right of all A^s using repeatedly (pTD , until it hits f2 an 
disappears (gives zero according to fl2"3p). This produces a sum of A" terms, and in each 
of them one of the fhj = 1 is changed to m^. Similarly we then remove p mN _ lt i, which 
turns in every term one of the remaining fhj = 1 to mjv-i, etc., until we remove p mi ,i 
and the remaining fhj = 1 in each term is turned to mi. Thus we used ( p7|) and ( f23"l) to 
remove all p's, which turns all = 1 to rrij, and we obtain AH different terms where we 
have the rrij in all possible different orders. The result after these manipulations can be 
written as 

$ = Ph,lh 2 ,2 ' ' ''%,iV 5 l,2,...,ivA,mi^2 1 m 2 • ■ ■ A^^Q 

where «Si t 2,...,iV means symmetrization over all indices rrij. In a similar manner we can 
remove all pe j, and this now turns the indices £j = j to tj. Since all ij are different and 
equal to only one of the corresponding indices of the A\ the number of terms does not 
increase and we obtain 

$ = S lt2 ,...,NA\ imi A\ 2in2 ■ ■ ■ A\ NmN VL = 
V A j A j ■■■A j Q 

PdS N 

which is the final formula for the eigenstate in this special case. 

It is easy to generalize this to the general case by inserting and into Eq. 



(pq). Again we can remove all the p's and p's using (|27|) and (P3|). Removing the p's 
turns all m^'s to m/s, but we obtain AJ x A2! x • ■ ■ x A^! terms corresponding to how 
often the numbers 1,2, ...,L appear in {fhj}^ =1 in (^4|): 1 appears Ai times and thus 
one has symmetrization in the first Ai indices rrij, 2 appears A2 times and thus one has 
symmetrization in the next A2 indices rrij, etc. Similarly, removing the p's turns the ^'s 
to £/s, but we further increase the number of terms by factors \\\ x \ 2 \ x ■ ■ ■ x \ N \ 
where A» is the number of times % appears in \lj\^ =x . The symmetrizations are now in the 
indices lj at the positions j where the ij coincide. The numbers Ai define another partition 
[A] = [Ai, A 2 , • • • , A^J of N, e.g. if [A] = [5, 3, 2] then [A] = [3, 3, 2, 1, 1] (readers familiar with 
Young tableaux will recognize this as the conjugate partition 0). The resulting formula 



for $ ( p6[ ) can be written conveniently in terms of partial symmetrization operators, 

<?• ■ ■ ... 4"!" O — 
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'Jp(l) '■32" l 3p<2) ^K" l Op(K) 

P&S K 



(this sum of K\ terms corresponds to symmetrization in the K indices rrij at the positions 
indicated by the subscripts), and another similar symmetrization but with respect to the 
indices 



9- • ■ 4* ... 4^ O •— 

c, JiJ2,-Jjf /1 « I m 1 /1 f 2 m 2 ^jvmjv ' — 

V A] m ■ ■ ■ A] ■■■A] ■ ■ ■ A\ ■ ■ ■ A) m Q. (39) 

P&S K 

The general formula is 

® = 5i,2,...,Ai'Sa 1 +1,Ai+2,...Ai+A2 " " " "Ai+...+Ajv-i+1,...,JV x 

xS 1 

,Ai+l,Ai +A 2 +l,...,Ai + ...+A Xi +l'52,Ai+2,...,Ai + ...+A^ 2 +2 • • • i j.^Ai+Aa+A-^ _ +Ai 

\N) (40) 

with \N) defined in (|T9]). To digest this general formula is is useful to consider few special 
cases, e.g. 

$ = «S 1)2i 3 J 4,5«56,7 J 8<S 9 ,io«5 1 , 6j9 «S 2! 7,io<S3,8|iV) for [A] = [5, 3, 2]. (41) 

(we used S4 = S$ = I) etc. A simple method to determine the subscripts in fl40"|) is to draw 
the Young tableaux corresponding to the partition [A] but write the numbers 1,2, ... ,N 
instead of the usual boxes, in increasing order from left to right and from up to down. 
For our example this gives 

1 2 3 4 5 

6 7 8 for [A] = [5,3,2]. 

9 10 

For each column write one S and for each row one S, and the subscripts of the S in (|40|) 
are given by the numbers in the columns and the subscripts of the S by the numbers in 
the rows. (The order in which the iS's or the <S's are written is irrelevant, of course.) 

To make contact with the results obtained in the next Section we note that the inter- 
action can be written also as sum of transpositions Ttjk) m another representation T of 
the permutation group Sn- Defining 

T(jk)A\ imi ■ ■ ■ A iNmN Q = A\ imi ■ ■ ■ A lkm , ■ ■ ■ A t . mh ■ ■ ■ A\ N>niN Vt (j < k) (42) 

(i.e. £j and £k are interchanged) we observe that T^^T^k) = —I (since this amounts to 
interchanging A\. m . and A\ m which is the same as multiplication with —1). We thus 
can also write 

H int \N) = -2g f mW)- (43) 

l<j<k<N 

We thus have two representations T and T of in our model which are, however, 
not independent but such that their product equals the representation P — > (— l)' p l (= 
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parity of P G Sjv). This fact allows us to somewhat simplify Eq. (|40|): we can replace 
the symmetrizers with respect to the indices £j by antisymmetrizers with respect to the 
indices rrij, i.e., the symmetrizers S are equal to 

<? ■ A^ A^ ... 4"!" O — 
~K~\ S ( _1 ) |P ' A li m i ' ' ' A lx^ p(1) ' ' ' A l 2 ™^ (2) ' ' ' A L m ^ ( x) ' ' ' A ^N,m N V. (44) 

Inserting that, we write $ = Y^\N). We now observe that is, up to normalization, 
the Young operator (or Young symmetrizer JFJ) associated with [A]: this result will be 
obtained in the next section using group theory. 



Remark: It is interesting to note that the model remains solvable if we add to the Hamil- 
tonian 7i = TCo + TCi nt Hartree-Fock type terms of the following form, 

H-HF = : {Vlrnpupmm + Vl m p~e£p mm + W tmPllpmm^ '■ (45) 

£,m 

where V, V and W are arbitrary model parameters: the eigenstates $ in Eq. ([HJ) remain 
the same, and the eigenvalues £ are changed by adding Yj&jWklj + ^» i m j + Wt imj ). We 
will come back to that in Section 9. 



6. Solution II. Group theory approach. As mentioned, there is a representation T 
of the permutation group Sn on the iV-particle states of the model, 

^ , \mi\m 2 ' ' ' A?jvmjv^ = A? imp{1) A? 2 m P{2) ' ' ' ^l N m P(N) ^ 

for all PeS N . 

As we now discuss, diagonalizing 7i amounts to decomposing T in irreducible repre- 
sentations (irreps). We first recall (p9|): the interaction applied to any state \N) defined 
in flT9|) is proportional to the sum of all transpositions, 

l<j<k<N 

to this state. The crucial fact allowing to solve the model is that the free Hamiltonian 7i 
commutes with all permutations Tp in ([46]) (since the eigenvalues So in fl28|) are invariant 
under all permutations of the indices rrij). If we therefore make an ansatz 

$ = apT P \N) (48) 

(with \N) short for any of the vectors in fliTS|)) and choose the real coefficients ap so that 

Cn apTp = c apTp (49) 

Pes N Pes N 
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for some real c, then $ is an eigenvector of 7i with eigenvalue 8 = 8 o + 2gc. 



Eq. (|49| ) can be interpreted as an eigenvalue equation for and can be solved as 
such. There are two obvious solutions: ap = 1 with c = N(N — l)/2 and ap = (— l)' p ' 
with c = —N(N — l)/2. They correspond the the eigenstates 

\ N h = E (±i) |p| t>W (50) 

Pes N 

with the corresponding energy eigenvalues 

£± = E ±N(N-l)g. (51) 

For fixed Sq these are the extreme eigenvalues. The total number of solutions is N\ 
generically. To appreciate the problem in Eq. (f49|) it is instructive to solve it by brute- 
force as discussed in Appendix B, but this approach is only possible for small N. 



To see that Eq. (f49|) is equivalent to a classical group theory problem, note that 
Cn = J2i<j is a class operator, i.e. it commutes with all permutations T P . Thus 
Cat is equal to a constant c in each irreps. Eq. (||9|) thus amounts to decomposing the 
representation T into irreps. The irreps of Sn are well-known. They are in one-to-one 
with the partitions [A] of N, and the value cpo of Cat in this irreps is well-known and equal 
to what we found in Eq. ( |3"3"D before (see e.g. Eq. (4-3) in [jl^jj ). If all Ej and all mj are 
different, then the multiplicity of these eigenvalue equals (&[a]) 2 where 

is the dimension the irreps [A] (this follows from the fact that T in this case is equivalent 



to the regular representation of Sn — see e.g. Theorem 3.25 and Eq. (4-4a) in [fL8|| ; the 
interested reader can find some examples in Appendix B). The corresponding Y2p a pTp 
is equal to the Young operator |6|, 0. We thus have recovered the result in the 
previous Section. The group theory argument is somewhat more powerful in that it also 
shows that no eigenfunction was missed. Moreover, the Young operators are projections 
up to normalizations: these normalizations are known and give the normalizations of the 
eigenstates. 



Remark: We thus found all eigenstates and eigenvalues of the model. One remaining 
problem is that we actually found too many: For a given \N), we obtained N\ eigenfunc- 
tions Y^\N), but all of them are linearly independent only if all iV ij& and all iV m/s 
are different. E.g. if all rrij = m are equal, then only one of these eigenvectors is non-zero 
(the others will all vanish due to the fermion anticommutator relations) and similarly for 
states where all Ej are equal. If we put identical E/s in groups then the ^-degeneracies 
can be characterized by the numbers Uj of elements in the different groups. This defines a 
partition [u\. For example, for (£j) = (2,5,5,2,3,2,5,2) one gets [u] = [4,3,1]. Similarly 
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the degeneracies of the m/s can be characterized by another partition \p]. Thus the de- 
generacies of the eigenvalues S depend on three partitions. To find all these multiplicities 
seems like a rather non-trivial exercise in combinatorics. If one knew these multiplicities 
m, one could compute the partition functions Z = Tracer exp(—/3H) of the model as 
follows, 

Z= Yl e-^w^+^m^je-^l (53) 

where (u) = (yx, 1/2, • • •) with vi > and [u] is the corresponding partition with the z/j's 
ordered (to have this sum finite one might still need a finite cut-off A, of course). 



7. Generalization to 2n dimensions. We now discuss the generalization to 2n dimen- 
sions. There one has (generalized) Landau eigenf unctions 0i )in which are labeled by 2n 
positive integers, 1 = (£x, ■ ■ ■ ,£ n ) and similarly for m. This is obvious in the coordinate 
system where the matrix B has Jordan normal, 



so that Hb in 
and X2j, j = l 
eigenfunctions - 



/ 





-Bx 



Bx 




V 





— B r 



\ 



B n 
J 



5„>0, 



(54) 



is the sum of n terms Hb depending only on the coordinates xj^x 
. ,n. Thus 0i, m ( x ) is just the product of n two dimensional Landau 
nj (x2j-i,X2j), and the corresponding eigenvalues of Hb are 



(55) 



Then all what we discussed for two dimensions goes through with £, m replaced by l,m 
etc. We now observe that we can map the vectors 1 = (£x, ■ ■ ■ ,£ n ) i n a one-to-one way 
to a single positive integer I (e.g. for n = 2 one such map is (1,1) — > 1, (1,2) — > 2, 
(2,1) -> 3, (1,3) -> 4, (2,2) -> 5, ...(£x,£-£x) -> £{£ - l)/2 + £x). Doing that, 
all formulas obtained in two dimensions hold true as they stand also in 2n dimensions 
with the only difference that the eigenvalues Eg and E m are given by somewhat more 
complicated expressions. However, our construction of eigenstates and eigenvalues above 
does not rely on the explicit form of these eigenvalues, and we thus have obtained all 
eigenstates and corresponding eigenvalues for arbitrary dimensions 2n. 



8. More general solvable models: An obvious generalization of our model would be 
to add a further interactions corresponding to a higher Casimir of gl^, for example the 
term H {3) = g (3) : Tr(p 3 ) :, i.e. 

: ^int : = 9 {3) ^2 ■ PkiPimPrnk ■ (56) 
k,£,m 
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with the colons indicating normal ordering as usual (i.e. move all A* to the right of all 
A's). Again we can write this interaction also in terms of the p's, tC^ 
This term comes the following 3-body interaction, 



: Tr(p 



n 



(3) 
int 



(3) 



(57) 



where regularization by normal ordering amounts to a renormalization of the 2-body 
interaction constant g and the chemical potential //. Applying this interaction to a state 
\N) (0) yields 

:^i ( n t ) :|iV)=^ (3) ^T ( , M) |iV) (58) 

where the sum on the r.h.s. is over all 3-cycles of the permutation group Sn 0. The 
eigenstates of the model which we constructed can be chooses such that they are also 
eigenstates of this interaction. The possible eigenvalues of this are known: they are 
characterized by partitions [A] and are given by (see e.g. Eq. (4-3) in 



A' 



.(3) 



1 ' 2 iV-^V 2 + $> 

1=1 



A 



3/ - - ) A, + - 1) 



One can can add further interactions TC^ oc : Tr(pP) 
solvable model. It is interesting to note that 



(59) 



V 



4, 5, . . ., and still have a 



: Tr(f) 



(60) 



9. Outlook and open questions. 



• More general models. The model solved in this paper is only the simplest in a class 
of similar models: one can add a flavor (or spin) index to the fermion operators and 
thus increase the symmetry from gl^ to g\ k ® gl^. This allows for additional types 
of interactions (spin-spin-like for k = 2, e.g.). Models of this kind can be obtained, 
e.g., by truncating the 2D Hubbard model: simplifying the 2D Hubbard interaction 
by keeping only particular terms |4j]. In 2D one can keep, in addition to Hartree- 
and Fock-terms (leading to mean field theory), also particular 'mixed' terms and 
still have a model which, as I believe, is exactly solvable H. 

• Mean field vs. correlations. To put my results in perspective, I recall a well-known 
class of exactly solvable models which can be defined by the Hamiltonian 7i = 

+ ^hf with the free part in (|i~6"D and the Hartree-Fock interaction in fl45D (this is 
only a special case 0). For this model all states \N) in fll9"D are eigenstates: Hartree- 
Fock interactions only change the energy eigenvalues but not the energy eigenstates. 
These eigenstates are Slater determinants: no correlations. For the models 7i = 
TCo + TCi n t + TCuf the eigenstates are highly non-trivial linear combinations of Slater 
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determinants (generically: there are exceptions, of course): for fixed E/s and m/s 
one has (generically) many different Slater determinants where the £/s and mj's are 
distributed over the fermions in all kinds of different ways, and only very particular 
linear combinations of these are energy eigenstates. For example, if all Ej's and 
all m/s are different, the 10-particle eigenstate in Eq. ( f4l|) is a particular sum of 
5!3!2!3!3!2! > 10 5 terms! I believe that most of these states cannot be written as 
Slater determinants (a proof of this would be very welcome). Thus, different from 
Hartree-Fock models, our model should describe correlated fermions. 

Phase transitions? In this paper I only demonstrate how to solve the models. To 
explore the physics of the solution is left to future work. However, to give a glimpse 
in that direction, I now discuss one special case of physical interest in which one 
can see by simple means that the solution can describe interesting physics: Assume 
n — 2, a — 1 and 6 = 0: 

U = ^ \B\ {m - l)pmm + 9^'- PimPmi ■ (61) 
m>l l,m 

(I set the chemical potential fi to some convenient value). As discussed, this is a 



toy model for a quantum Hall system |17 |. For b = it is natural to keep the 



cutoff A finite: it has a natural interpretation as a spatial cutoff, and v = A/A is 
the fermion density ('filling factor') of the QH system. The energy eigenvalues are 
sums of two terms: the kinetic energy £q = \B\ Ylf=i( m j ~ 1) an d the correlation 
energy £ corr = 2gc[x\. Let N < v. Then obviously £ o = is minimal if all rrij = 1 
(i.e. all fermions are in the lowest Landau level). Then all £j need to be different, 
and (for fixed £/s) there is a unique eigenstate \N) with correlation energy £ COTT = 
gN(N — 1). This obviously leads to a minimal energy £ = gN(N — 1) if g < 0, but 
for positive g it can be preferable have some fermions with mj > 1: this increases 
the kinetic energy but allows to decrease the correlation energy. For example, if 
rrij = j (all different: we put each electron in another Landau level) and Ej = i (all 
the same, to have a simple specific example), \N) is eigenstate with total energy 
£' = N(N- 1)| B\/2-gN(N-l). Obviously, £' will be lower than £ for sufficiently 
large g. Thus, if one increases the coupling from g < zero to g > \B\/A there 
must be some drastic change of the ground state in between — possibly a phase 
transition? This QH model becomes more interesting if one adds Hartree-Fock terms 
as in Eq. (f45|) : this allows to lift the degeneracy (fermions in the same Landau level 
can repel each other, e.g.). 

Meaning of dimension? From an abstract point of view, the model we solved looks 
the same in all dimensions 2n+ 1. This might seem somewhat puzzling, in particular 
if one recalls that renormalizability of quantum field theory models usually very 
much depend on dimensions. This does not seem to be the case here: these models 
therefore seem to challenge our usual expectations (and prejudices) about locality 
and dimension in quantum field theory. I should add that this is less unusual from 
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a solid-state point of view: the above-mentioned Hartree-Fock models also look the 
same in all dimensions. 
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Appendix A. The star product of the Landau eigenfunctions. Fact 1 and Fact 2 
in the main text are known since a long time. However, due to the weight they carry in 
my discussion I felt I should also include an (elementary) proof. This is the purpose of 
this Appendix. 

We assume 2n = 2 and write H B = H i = Pf + P 2 2 and H-b = H 2 = P 2 + P 2 where 
Pi = —id\ — \B\x2, P2 = —id 2 + \B\xi 

p x = -id x + 151x2, P 2 = -id 2 - \B\xt (62) 
with \B\ > 0, dfj, = q|jt and (x M ) = (x 1 , x 2 ) coordinates on R 2 . We observe that 



give a representation of the Heisenberg algebra [p M , q v ] = —i5^ u etc. Thus H% = 2\B\(pl + 
q\) is just a harmonic oscillator Hamiltonian which (on the Hilbert space of functions 
in two variables) is highly degenerate. The operator H 2 = -Pj 2 + P 2 = 2\B\(p?, + q%) 
is the 'complimentary harmonic oscillator'allowing us to resolve this degeneracy. The 
Landau eigenfunctions are the common eigenfunctions of H\ and H 2 . To construct them 
we introduce creation-and annihilation operators = -^(TWn + %) obeying the usual 
commutation relations. Then = 4\B\ (a+a~ + |). The common eigenfunctions of H 
and H 2 are therefore^ 



where |0) = and £, m positive integers. The eigenvalues are 4|P|(£— ^) and 4|P|(m — 
|), respectively. 

5 The phase factor (— i) m_1 is inserted for convenience. I find it also convenient to label these states 
by positive integers (i.e. what I call t — 1 is usually called £). 
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We now compute the normalized Landau eigenfunctions in position space, 

<^ m (x) = (x\£,m) (65) 
where x = (x 1 ,^ 2 ). It is convenient to define 

z = Xi + ix 2 , z = Xi — ix 2 (66) 
and to introduce the generating function 



F a , t (x) = -7^=f^+i,m+i(x) = 
o n v v-m\ 

£,m=0 

= E -^^<^l(«i) £ («2) m |0> = (z,z|e^-^|0). (67) 

^,m=0 

Since z = (a* + ia 2 )/-\/TB[ and z = (ai — ia^/VT^I one S e ^ s 

_p st (x) = ^ x |e sa i _i * a 2e tai+,isa2 |0) = (x|eV^I( sz+t ^e' sa i~ i * a 2' iai+,i ' sa2 ] //2 |0) = 

= ev^(-+ t5 ) e - s *0i,i(x). 

The normalized ground state wave function can be computed by solving a^i^ = 
with ai = (9i — ?<9 2 + | -B |(^i — ix2))/2y/\B\ and similarly for a 2 . One thus finds 



/■„(x)= ,/M e -^ eV ^(-+^)-l B H-l 2 /2. (68) 



Using this generating function it is easy to check that the Landau eigenfunctions thus 
defined are a complete orthonormal basis. 

To prove Fact 1 we compute F Sljtl * F S2jt2 assuming 9 = B^ 1 . We use 

(f*g)(x) = (2tt)- 2 / d 2 k ! d\f(k)g(q)e i W~ 1 ^ q2 - k2 ^e i ( k+ *> x (69) 

where / is the Fourier transform of /. We first compute (K — k± + ik 2 ) 



F s t (k) = [ — e- ik ' x F 8t (x) = J -!—e st e HsK + tR)/y/W-\K\V2\B\ 
(we computed a Gaussian integral). Thus 



F si , tl ★F S2it2 (x) 



(27r)- 2 / rf 2 k / ^2 qe -i0(fc ig2 -fc 292 )-i(k+q).x 

Jr 2 Jr 2 



1 s 1 t 1 +s 2 t 2 J(s 1 K+t 1 K+S2Q+t2Q)/J\B\~(\K\ 2 +\Q\ 2 )/2\B\ 



x— — e 1 1 1 ^ " e 
tt LB 



27 



e -sit 2 +s 2 ti e (siz+t 2 2) v /|B|-|B||z| 2 /2 
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(again a Gaussian integral), i.e. 



47T 

equivalent to 



F Butl *F aatta = W^e^F^x) (70) 



i,mi *4>e2,m 2 — \J \^~~;fimij2<f>£i,ms> ■ (^1) 

and completing our proof of Fact 1, including the normalization. Our discussion in the 
beginning of this Appendix also provides a proof of Fact 2. 



Appendix B. Pedestrian solution of Eq. (|49|). In the main text we give the general 
result for all N. To appreciate this result and have a few special cases, I present here a 
brute-force solution for small N. 

One can represent Cn (W?h by a N\ x AM-matrix in the following way. Define a 
Hilbert space isomorphic to ~R Nl by introducing the orthonormal basis \P), P E Sn- 
Then 7Tp(P', P") := (P', P ■ P) defines a representation P — > Tip by iV! x iV! matrices (the 
reader familiar with group theory will recognize that it is just the regular representation 
of 7r). We thus can solve (|49|) by constructing and diagonalizing the AH x A!-matrix 
Cn := $^i<3<k<jv 7r (j' fe ) - We now describe the results for iV < 4 and give a few examples 
for eigenfunctions. 

For N = and iV = 1 it is easy to find all eigenstates and eigenvalues of Ti,, 

HVI = (N = 0) (72) 

and 

HA\jl = (E e + E m )A\jl (N = 1). (73) 

The first non-trivial case is iV = 2. Labeling the basis \P), P G S%, in the following way, 
1 12) = 1 and |21) = 2, (we use an obvious short-hand notation for permutations on the 
l.h.s. of these equations) it is easy to see that 

*-(;;)• 

which has eigenvectors (aj, a§) = (1? il) corresponding to the eigenvalues 

c=±l (N — 2). (74) 
From that one obtains the following eigenstates of 7Y, 

|2) ± = (a t I + a 2 T (12) )|2) = (< mi 4 2m2 ± AmA^) « (75) 

(we identified Tj = T 12 = J and = T 2i = T( 12 )) corresponding to the eigenvalues 
£ = £ ± 2g. Note that if m 1 = m 2 then \N) + = and if l x = £ 2 then |JV)_ = (due 
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to the Pauli principle). Thus there are two independent eigenstates only if mi 7^ m-i and 
h 7^ otherwise there is only one. Still, the N = 2-eigenstates thus constructed provide 
a complete orthonormal basis in the N = 2-subspace of our fermion Hilbert space. 

For iV = 3 we label the basis as |123) = i, |312) = 2, |231) = 3, |213) = 4, 1 132) = 5, 
and 1 32 1 ) = 6, and we obtain 



1 1 
1 1 
1 1 

111 
111 
111 



(the matrix elements not written are zero). The eigenvalues of this matrix are 

c= 3(1), 0(4), -3(1) (N — 3) (76) 

with the numbers in the parenthesis indicating the multiplicities. It is not difficult to also 
write down the corresponding eigenstates. For example, (aj, . . . , ag) = (1, —1, 0, 0, 0, 0) is 
an eigenvector of C3 with eigenvalue c = 0, and the corresponding eigenstate of TC with 
eigenvalue £ = £ is 

(T 12 3 - I312) |3) = [Al imi A\ 2m2 A\ 3m3 ± A\ im3 Al 2mi A\ 3m2 j Q (77) 

etc. If all three m^s are different and also the three i^s, these eigenstates all are linearly 
independent, but otherwise the number of linearly eigenvectors can be less than 6 and not 
all eigenvalues £o + 2gc are realized. I also constructed and diagonalized the 24 x 24-matrix 
C4 and found the following eigenvalues (and multiplicities), 

c = 6(1), 2(9), 0(4), -2(9), -6(1) (N = 4). (78) 

Again, the corresponding 24 eigenstates are linearly independent only if the four rn^'s are 
all different and the same for the IjS. At this point this brute-force approach clearly 
becomes too cumbersome. Fortunately for us, the eigenvalues and eigenvectors of all the 
matrices Cn are all known from group theory. 
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